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We calculate the diameter and chirality dependences of the binding energies, sizes, and bright-dark 
splittings of excitons in semiconducting single- wall carbon nanotubes (SWNTs). Using results and 
insights from ab initio calculations, we employ a symmetry-based, variational method based on the 
effective-mass and envelope-function approximations using tight-binding wavefunctions. Binding 
energies and spatial extents show a leading dependence with diameter as l/d and d, respectively, 
with chirality corrections providing a spread of roughly 20% with a strong family behavior. Bright- 
dark exciton splittings show a l/d 2 leading dependence. We provide analytical expressions for the 
binding energies, sizes, and splittings that should be useful to guide future experiments. 
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Diameter and chirality trends are one of the most use- 
ful concepts in nanotube science. Often, new physics 
arises when the diameter and chirality dependences of a 
given property are fully disclosed. A classic example is 
the analysis of "family patterns" in optical transitions 
combined with the diameter dependence of vibrational 
frequencies that paved the way to reliable (n, m) assign- 
ments of single-wall carbon nanotubes (SWNTs) [lj and 
posed the fundamental "ratio problem" j2j. Therefore, 
a reliable determination of diameter and chirality trends 
of a given nanotube property, even when this is accom- 
plished by simplified models, is often as important as de- 
termining accurately that property for a limited number 
of tubes. Moreover, when a reliable model for trends is 
coupled with an accurate ab initio theory that determines 
its parameters, the model acquires quantitative and pre- 
dictive powers. 

The exciton concept solved the "ratio problem" 0, 
and it is now widely accepted that the optical spectra 
of carbon nanotubes are dominated by exciton features 
[3 IE El • Recent experiments based on two-photon 
spectroscopy E3 and Raman spectroscopy on elec- 
trochemically doped samples have provided the first 
experimental evaluations of exciton binding energies for a 
few single- wall carbon nanotubes (SWNTs). However, a 
full description of diameter and chirality dependences of 
exciton properties in SWNTs has not yet been provided, 
neither experimentally nor theoretically. Ab initio calcu- 
lations are restricted to a few small-diameter tubes 00. 
Perebeinos et al. |fj have extracted scaling relations of 
binding energies and sizes with diameter from model cal- 
culations, but the chirality dependence has not been ad- 
dressed. Semi-empirical calculations have also been done 
for a larger variety of tubes 0, but again systematic 



diameter and chirality trends have not been extracted. 
Finally, the important issue of bright-dark exciton split- 
tings have been addressed in detail by considerably fewer 
calculations [HElEi- 

In this work, we calculate the full diameter and chi- 
rality dependences of exciton properties in SWNTs. 
We employ a symmetry-based, variational, tight-binding 
method, based on the effective mass and envelope func- 
tion approximations |lfj Hfij | . Since we explictly impose 
the symmetry of the exciton wavefunction, we can calcu- 
late properties of bright and dark excitons. Our model is 
parametrized by ab initio results. We calculate binding 
energies and sizes for the lowest-energy bright excitons 
(those usually associated with the En singularity in the 
single particle joint density of states), as well as dark- 
bright exciton splittings for a large number of SWNTs. 
From these results, we extract reliable analytical expres- 
sions for the diameter and chirality dependences of such 
properties. 

Our variational exciton wavefunction is written as: 



(1) 



where c (r e ) and (j> v (rh) are conduction (electron) and 
valence (hole) single-particle states. The sum is restricted 
to the four band-edge states c = ±m and v = ±m 
from the top of the valence and bottom of the con- 
duction bands. Note that both valence and conduction 
band edges are 2-fold degenerate for both zigzag and 
chiral tubes, once time-reversal symmetry is considered 
[Til Hsf . The single-particle wavefunctions arc labeled 
by their quasi-angular momentum quantum numbers +m 
and — m and they are taken from properly symmetrized 
wavefunctions of graphene expanded in a 7r-orbital tight- 
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TABLE I: Symmetries, degeneracies, optical activities and 
coefficients A vc for excitons in zigzag and chiral tubes. The 
symmetries are described by the irreducible representations 
in both group of the wave vector 0] and line group |l7|(in 
parenthesis) notations. The label m is the quasi-angular mo- 
mentum quantum number of the double-degenerate exciton. 
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three terms: direct, exchange and kinetic energies. Here, 
we treat singlet excitons only. The direct term is written 
as: 

<K d >= f r(re,f h )VS cr (fe~r h )ij(f e ,r h )dr e dr h 



vc,v'c' §.i,R 2 

X e- l -^^U^ no {\Ri-R2\)m 

where Vq ct is the screened Coulomb interaction and we 
wrote the direct energy in terms of the tight-binding ex- 
pansion coefficients of the single-particle wavefunctions 
in a p z -orbital basis (p(f — Ri) centered in the atomic 
positions Rf. 



= Ycn{Rj)f(r- Ri)- 



(3) 




z, - z, (nm) 



FIG. 1: (Color online) Lowest-energy singlet exciton wave- 
function of the (11,0) tube. Black thin line: Ab initio \i/j\ 2 
after integrating out on the coordinates perpendicular to the 
tube. Red thick dashed line: Envelope fit using a gaussian. 
Blue thick line: Envelope ht using a Whittaker function. No- 
tice that the two fits are almost indistinguishable. 



binding basis. The coefficients A vc , responsible for 
the quantum interference between pair excitations, are 
then completely determined by symmetry, as described 
in Table [Q 

We choose a gaussian envelope function. This choice 
is justified both by a fit of the ab initio exciton wave- 
functions 0, as shown in Fig. ^ and by an analogy 
with the regularized Coulomb potential problem in ID 
for which the ground-state Whittaker function 
closely resembles a gaussian. The gaussian width or ex- 
citon size a is the only variational parameter in the prob- 
lem. The constant C normalizes the exciton wavefunc- 
tion: / J \ip(r ei r h )\ 2 dr e dr h = 1. 

We minimize the exciton energy that is composed of 



The Coulomb integrals between sites are parametrized 
by the Ohno formula |2*l) : 



Uohno(R) 



U 



^(^u Ry + 



(4) 
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The onsite Coulomb repulsion C/ = 16 eV and the di- 
electric constant e = 1.846 are chosen to reproduce the 
ab initio values for the binding energy and bright-dark 
exciton splittings for the (11,0) tube and kept constant 
for all other tubes. The exchange energy is given by: 



< K x >= 2 J ip*(f e ,r e )Vc(re - r h )^[f h ,r h )dr e drh 
= 2C 2 Y A *vc A v>c Y c v {Ri)c c (Ri)c:,(R 2 )c* c ,(R 2 ) x 

vc,v>c' Ri,R 2 

xU hno(\Ri-R2\)-{5) 

In this case, the unscreened Coulomb interaction Vc is 
parametrized by taking e = 1 in Eq.QJ. Finally, the 
kinetic energy associated with the exciton relative coor- 
dinate is simply that of a gaussian envelope: 



< T >= 



Am*a 2 



(6) 



where the exciton reduced mass m* is given by 1/m* = 
l/m e + 1/m/j. We use the diameter- and chirality- 
dependent electron (m e ) and hole (rrih) effective masses 
obtained from tight-binding calculations [22J. 

To test our model, we compare in Table [H] our varia- 
tional binding energies with ab initio ones obtained from 
solving the Bethe-Salpeter equation 3] for a few zigzag 
tubes. The agreement is excellent, except for the E22 ex- 
citon in the (7,0) SWNT. This discrepancy can be under- 
stood: For the small-diameter (7,0) tube, the E22 exciton 
size becomes extremely small (a = 5.3 A) and therefore 
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TABLE II: Ab initio and model binding energies for bright 
En and E22 excitons for a few small-diameter SWNTs. 





E h 


(eV) 


Bf (eV) 


Tube 


Ab Initio 


Model 


Ab Initio Model 


(7,0) 


0.89 


0.87 


1.13 1.61 


(8,0) 


0.99 


1.03 


0.86 0.92 


(10,0) 


0.76 


0.68 


0.95 1.09 


(11,0) 


0.76 


0.76 (fitted) 


0.72 0.75 
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FIG. 2: Binding energies for the lowest-energy bright excitons 
in 38 SWNTs with varying diameter and chirality. The dots 
are our model results and the red lines represent the analytical 
fit using Eq.JTJ. The labels indicate the (2n + m) families. 



the envelope-function approximation is not expected to 
be valid in this regime. Notice also that our model cor- 
rectly captures the family oscillations in the binding en- 
ergy. 

Fig. |2Jshows the binding energy of the En bright exci- 
ton as a function of diameter for 38 SWNTs covering the 
full range of chiralities. The (2n + m) family indices are 
indicated in the figure. One clearly sees the well-known 
family pattern reminiscent from the so-called Kataura 
plots for optical transition energies, As expected, 

binding energies decrease with tube diameter. Chirality 
effects are also strong, contributing to about 20% spread 
in the binding energies for a given diameter. It is clear 
that excitons in (2n + m) mod 3 = 1 (MODI) tubes have 
generally larger binding energies than in (2n + m) mod 
3 = 2 (MOD2) tubes. 

Fig. |3|shows the exciton sizes as a function of diameter. 
Again, as expected, exciton sizes increase with diameter 
and they show the opposite MOD1-MOD2 trends as com- 
pared to binding energies. Notice that even for tubes as 
small as 0.5 nm in diameter the En exciton sizes are al- 
ready several times larger than the carbon-carbon bond, 
thus justifying the use of the envelope-function approxi- 
mation. 
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FIG. 3: Sizes of the bright lowest-energy excitons for 38 
SWNTs with varying diameter and chirality. The dots are 
our model results and the red lines represent the analytical 
fit using Eq.JTJ. The labels indicate the (2n + m) families. 



Analytical expressions for diameter and chirality de- 
pendences, although sometimes lacking a deeper physical 
justification, can be extremely useful for a quick evalua- 
tion of a variety of nanotube properties. We succeeded 
in finding simple yet very accurate analytical approxima- 
tions for both binding energies and sizes: 



1 / B 
E b = -[ A+- + C£,+D£ 2 



d\ d 

a = d(E + FZ + Gf), 



(7) 



where d is the tube diameter in nm and £ = {—l) u cosiO / d 
captures the chirality dependence jl^. The best fits are 
given by A = 0.6724 eV.nm, B = -4.910 x 10~ 2 eV.nm 2 , 
C = 4.577 x 10~ 2 eV.nm 2 , D = -8.325 x 10~ 3 eV.nm 3 , 
E = 1.769, F = -2.490 x 10" 1 nm and G = 9.130 x 10~ 2 
nm 2 . These analytical fits are plotted in red lines in 
Figs. |2 and |3 together with the numerical results. The 
agreement is nearly perfect. 

Our theory also allows for an estimation of chirality 
and diameter dependences of exciton splittings among 
exciton states of the ground-state complex of the same 
Eii. These splittings are fundamental to understand a va- 
riety of optical properties of carbon nanotubes, such as 
the quantum efficiency for light emission and the exciton 
radiative lifetime (l2T.ll3|. We find that the lowest-energy 
exciton for all SWNTs is the singly- degenerate dark state, 
due to its vanishing exchange energy [12J. Defining the 
exciton splittings from the lowest-energy exciton to the 
bright exciton and to the double-degenerate dark exci- 
ton as Si and 82, we find the following dependence on 
diameter and chirality: 

6i = ^ (Ai + + Cde) , (8) 
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with A x = 18.425 meV.nm 2 , B x = 12.481 meV.nm 3 , 
Cx = -0.715 meV.nm 3 , A 2 = 32.332 meV.nm 2 , B 2 = 
7.465 meV.nm 3 and C 2 — —2.576 meV.nm 3 . So, in dis- 
agreement with Perebeinos et al. [l3l | , we find the leading 
dependence of bright-dark splittings to be 1/d 2 . This is 
precisely the dependence of the exchange energy < K x > 
on diameter. 

It is instructive to explain on physical grounds the lead- 
ing dependences on diameter of the exciton sizes, bind- 
ing energies and the bright-dark splittings. The exciton 
sizes a scale like d because the ID Coulomb potential is 
smoothed out or regularized over the scale of the tube 
diameter d and this sets the length scale of the bound 
state (recall that in a pure ID system with no lateral 
size, the Coulomb potential gives a delta-function ground 
state with infinite binding energy). The binding ener- 
gies go like 1/d because a scales like d and Coulomb 
interactions go like inverse distance Q. The scaling of 
dark-bright splittings mirrors the scaling of the exchange 
energy < K x > which goes like 1/d 2 because < K x > is 
the self-interaction of a neutral charge distribution with 
dipole moments: The long-range part (from distances 
larger than d) can be written as J d dx/x 3 ~ 1/d 2 . 

We now compare our results to the available experi- 
mental determinations of the exciton binding energies to 
date. Two-photon spectroscopy have been performed for 
SWNTs in a polymeric matrix B, 0] and in D2O solu- 
tion wrapped by a surfactant [lQj. These environments 
should provide extra screening, so these results should 
not be directly compared with ab initio theory for iso- 
lated tubes. However, in our variational scheme, it is 
very easy to investigate the influence of screening and 
to adjust the dielectric constant e to match the exper- 
imental results. In fact, we find that binding energies 
follow very nicely the scaling E\, oc e~ 1A proposed by 
Perebeinos et al. j|. Therefore, it is straightforward 
to apply Eq.J7|) for SWNTs in any environment, pro- 
vided that one scales the binding energies by using the 
appropriate phenomenological dielectric constant. For 
instance, taking e = 3.049 gives binding energies in ex- 
cellent agreement (standard deviation of 0.02 eV) for all 
13 SWNTs measured by Dukovic et al. Q. Similarly, the 
results of Maultzsch et al. [hJ for 6 different SWNTs are 
reproduced with a standard deviation of 0.03 eV using a 
slightly larger dielectric constant e = 3.208. 

In another recent experiment, Raman spectroscopy un- 
der electrochemical doping was used in nanotubes coated 
with a surfactant to give 0.62 eV and 0.49 eV for the 
binding energies of excitons associated with E22 transi- 
tions in the (7,5) and (10,3) SWNTs, respectively. [ll| 
We have also calculated the binding energies of those 
excitons. It should be noted that for E22 excitons, the 
MOD1-MOD2 oscillations in the binding energies are in- 
verted, i.e., MOD2 tubes have larger binding energies 
than MODI tubes of similar diameter. In fact, in dis- 
crepancy with experiment, we find that (7,5) and (10,3) 



tubes should have E22 excitons with similar binding en- 
ergies, even though the latter has a larger diameter. By 
using e = 2.559 we find the best possible "average" agree- 
ment with experiment: 0.54 eV for the (7,5) and 0.55 eV 
for the (10,3) nanotube. 

In conclusion, we have determined the full diameter 
and chirality dependence of exciton binding energies, 
sizes, and splittings in semiconducting SWNTs. All these 
exciton properties have strong diameter and chirality de- 
pendences, with a distinct family behavior. Comparisons 
between theoretical and experimental binding energies 
should be exercised with care, by ackowledging environ- 
mental screening effects. Our results should provide an 
useful guide to the interpretation of recent and future 
experimental determinations of exciton binding energies 
and other properties. 
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